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APPLICATION OF A TRANSONIC SINILARITT RULE TO 
CORRECT THE EFFECTS OF SIDEUALL BOUNDARY LAYERS IN 
TWO-DIMENSIONAL TRANSCNIIC WIND TUIWELS 
by 

William G. Sewall 
ABSTRACT 

A transonic similarity rule uhlch accounts for the effects of 
attached slde'.fall boundary layers is presented and evaluated by com- 
parison with the characteristics of airfoils tested in a two-diaienslooal 
transonic tunnel with different sidewall boundary-layer thicknesses. 

The rule appears valid provided the sidewall boundary layer both mains 
attached in the vicinity of the model and occupies a small enough 
fraction of the tunnel width to preserve sufficient two-dimensionality 


in the tunnel. 
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CHAPTER 1 


INTRODUCTION 

Since the development of wind tunnels, extensive attention 
has been devoted to the interference on wind-tunnel models caused 
by the tunnel wall boundaries. This interference is caused by 
the alteration of the streamlines near the wall from their free- 
air positions. Wind-tunnel interference in both two- and three- 
dimensional facilities has been addressed with analysis and 
facility modifications to reduce or eliminate it. 

In the past, the primary interest in two-dimensional wall 
interference concerned the upper and lower wall effects. Several 
of the linear analytical methods that have been developed to 
account for these effects are presented in reference 1. Some 
nonlinear methods have also been developed for transonic two- 
dimensional tunnels and are desfibed in references 2 and 3. 
Attempts to reduce these interference effects have resulted in 
facility modifications by making the upper and lower walls with 
either longitudinal slots, porous surfaces, or adjustable contours. 

The interference effects caused by the sidewalls in two- 
dime asional wind tunnels occur because of the presence of the 
sidewall boundary layers. The interference of attached sidewall 
boundary layers in two-dimensional tunnels results in a modifica- 
tion of the continuity equation because the geometric tunnel width 
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is effectively reduced by twice the sidewall boundary-layer dis- 
placement thickness. The two sidewall interference problems which 
have received the most attention are the growth of the sidewall 
boundary layer due to the shearing stress at the sidewall and the 
separation of the sidewall boundary layer due to interaction with 
large model-induced pressure gradients. The problem of boundary- 
layer growth due to shearing stress is accounted for in some wind 
tunnels by a slight outward inclination of the walls, and the 
problem of the sidewall boundary-layer separation can be con- 
controlled to some extent with suction or tangential blowing on 
the sidewall. 

This study concerns the Intermediate problem of the attached 
sidewall boundary-layer interaction with the pressure field of 
the model at transonic speeds. Earlier methods of accounting for 
this effect have been proposed for incompressible flow, as 
described in references 4 and 5. These methods considered the 
effect of the sidewall boundary layer as a change in the circula- 
tion about the airfoil. 

for subsonic and tiansonlc compressible flow, this effect 
can be formulated into similarity rules of the sidewall boundary 
layer to the model-induced pressure field. The analysis presented 
in this study applies elements of the derivation of the similarity 
rule given in reference 6 to the von Karman transonic similarity 
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rule. Experunental results from three airfoil tests, each conducted 
with varying sidewall boundary- layer thicknesses, are also presented. 
These results are used to evaluate the validity of the similarity 
rule at transonic speeds. 



1 


4 


CHAPTER 2 

EXPERIMEHTAL APPARATUS 
Facility and Test Conditions 

The effects of tho sidetnll boundary layers in subsonic and tran- 
sonic two-dimensional tunnels were investigated in the langley 6- by 
19-Inch Transonic Tunnel, presented in figure 1. This facility, 
described in detail in reference 7, is essentially a blowdown tunnel 
that operates at Mach numbers ranging from 0.3 to 1.0 with corresponding 
unit Reynolds numbers of S.O million to 7.S million per foot. 

The tunnel axis is oriented vertically with the flow direction 
upward, as shown in figure 1. The test section has solid, parallel 
sidewalls and slotted walls Joining the sidewalls that minimize the top 
and bottom wall interference mentioned in reference 1. 

This tunnel is configured for testing airfoil models, whicn span 
the tunnel, as shown in figuie 1, and have constant cross sections con- 
sisting of the airfoil shapes. Each end ol the model mounts into a 
turntable that fits flush into the sidewall of the test section. The 
turntables rotate together, allowing citanges in the model angle of 
attack. A movable rake mounted behind the model is equipped with four 
total-head tube probes and traverses the wake of the model to obtain 
wake total pressure measurements used in the drag calculation method 
from reference 8, 

Static pressure measurements both in the test section and in the 
contraction region 45.7 cm upstream of the start of the test section 
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are provided by a centerline row of orificea in the surface of one side- 
wall. On the same sidewall, a turntable with 17 orifices is used to 
obtain the Mach number at the model station when running tunnel empty 
for tunnel calibration. 


Models 

A photograph of two typical models tested in the facility is shown 
in figure 2. Both luxiels are constructed of stainless steel, have a span 
of 15.72 cm, and are Instrumented for pressure distribution tests. The 
15.72-cm chord model, which was the type used in this experiment, has 
rectangular tangs machined on the ends of the model to transfer the 
model aerodynamic loads to the model support system. The tubes have 
been placed inside the model and the coverplate has been welded in place 
over the tubes. The orifices are located in chordwise rows near the 
midspoii of the model and have a diameter of 0.33 mm. The model contour 
accuracies for these experiments were within ±0.013 mm. 

Three airfoil shapes were used for the experiment. The first model 
was a NACA 0012 airfoil model, which is a symmetrical airfoil that had 
been tt<».ed in many transonic facilities. The second airfoil to be 
I tvd was the supercritical SC-27 airfoil that represented the modern 
class of transonic airfoils. This particular airfoil has been exten- 
sively tested in two neighboring facilities and gave the first indica- 
tion of an interference problem that was thought to be caused by the 
sidewall boundary layer. The last model tested was that of the NLR-1 
airfoil, another supercritical airfoil that was designed for rotorcraft 
applications and wa*^ tcsiod in a neighboring facility. 
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All of Che airfoil models were tested trith a spanwlse strip of 
carburundum grit oa both upper and lower surfaces to force the laminar-' 
to-Curbulent boundary- layer cransicion. For the NACA 0012 airfoil, the 
transition particles were located at 0.075 x/c and were 0.089 nm In 
nominal height. The SC-27 and NLR-1 airfoil models both had Che transi- 
tion particles located at 0.05 x/c with a ntmiinal particle height of 
0.074 ram. The particle sizes were derermlned by reference 9 so as to 
insure sufficient size for complete transition without contributing 
extra protuberance drag. 

Artificial Thickening of the Sidewall Boundary Layer 
The sidewall boundary effects on two-dimensional airfoil testing 
were studied by examining results of several airfoil tests conducted 
with successively Chickened sidewall boundary layers. The sidewall 
boundary layers were thickened with a device similar Co that investi- 
gated in reference 10, which consisted of thin plates, each having 

three rows of pins protruding from Che surface, as shown in figure 3. 

One pl.ite was mounted on each sidewall in the tunnel contraction region 

(fig. 1), at a station 114 cm upstream of Che model leading edge as 

indicated in figure 4. The thickness of the sidewall boundary layer 
was controlled by the distance that the pins protruded from the plate 
surface. Three pairs of plates were used: the first pair had no pins, 

the second pair had pins extending out 2.54 cm, and the third pair had 
pins extending out 3.80 cm. The plates were mounted as far upstream as 
was conveniently possible to allow the wakes from Che individual pins 
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to be sufficiently mixed so that no peculiar behavior would exist in the 
sidewall boundary layers close to the model station. 

The sidewall boundary layer was surveyed along the test section 
centerline using total-head fixed rake-tube probes. These probes con- 
sisted of two rows of 0.76 am o.d. tubes and are shown in figure 5. The 
tubes %fere positioned to survey out to S.IO cm from the sidewall surface. 
The tube closest to the sldemll rested on the sidewall and was used to 
determine the skin-friction coefficient from the Preston tube calibra- 
tion in reference 11. The static pressure at the rake-tube probe posi- 
tion on the sidewall was determined from the static pressure measured at 
the sidewall location nearest to the probe location without the probe 
inserted. This static pressure was calibrated against the freestream 
Mach number so that with the probe mounted, the boundary-layer static 
pressure was determined from the freestream Mach number. 

The velocity ratio, U/Ug, iras calculated at each tube position on 
the rake-tube probe for the different sidewall boundary-layer thick- 
nesses. First, the calibrated static pressure was assumed constant 
through the boundary layer, and with the local total pressure at a 
particular cube, the local Mach number was calculated. Next, a static 
temperature distribution through the boundary layer, given by the 
equation 



( 1 ) 
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frcMii references 12 and 13, was used to calculate the local speed of 
sound. The local velocity at each tube location uas then determined 
from both the local Mach number and speed of sound. 

After the local velocities through the boundary layer were obtained, 
the boundary- layer thickness was determined from a least-squares regres- 
sion using the power law 



With the values of 6 and N and the static temperature given by 
equation (1), integrations were performed to get the displacement thick- 
ness, 5*, and the momentum thickness, 6. Details of these procedures 
are given in Appendix A. 

To examine the similarity of the thickened sidewall boundary 
layers, the velocity profiles at the model station were compared. 

Figure 6 shows the results of this comparison where the variation of the 
local velocity ratio, U/Ug, with the nondimens ional height, z/6*, 
.ippear for all three boundary- layer thickening configurations. Th,_je 
profiles were considered to match quite well. 

Another characteristic investigated for the thickened boundary 
layers concerned the relationships between the velocities in the 
boundary layer and the wall shear stress. One relationship for turbu- 
lent boundary layers in a zero pressure gradient has been established 
with the law of the wake of reference 14. The correlation of the three 
artificially-thickened boundary layers with the law of the wake is shown 
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in figure 7. The c''rrelatlon is reasonable except for the inner regions 
of the boundary layers. This is as expected because the law of the wake 
applies to the outer regions of the boundary layer while the more 
commonly known law of the wall applies to the inner region. The rake- 
tube boundary- layer probes used appear to define primarily the outer 
region of the boundary layers. 

From the law of the tiake correlation in figure 7 and the similarity 
of the velocity distributions in figure 6, it is concluded that the 
three boundary-layer Chlckenlrig configurations produced turbulent 
boundary layers similar to those on a smooth fla plate. 

The artificially-thickened sidewall boundary layers were also sur- 
veyed In the vicinity of the model station, tunnel empty, to find the 
variation of boundary- layer displacement thickness. Figure 8 shows this 
variation to be small for all three boundary-layer thickening configura- 
tions. This slow growth rate of the sidewall boundary layer is similar 
to that due to the shearing stress on a long, flat plate. 

The value of the shape factor, H, ranged from 1.39 to 1.59 for 
the three boundary- layer thickening configurations. 
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CHAPTER 3 

ANALYSIS 
Primary Coacepts 

Consider steady, isentropic, small-perturbation flow In a nominally 
two-dimensional airfoil wind tunnel. Let the Cartesian coordinates In 
the freestream, normal, and spanwlse directions be x, y, and s; and 
the respective velocity components be U, v, and w, as shown in fig- 
ure 9. The effective tunnel width is b-25* where b and 6'* can 
vary slightly with respect to x and y, and the boundary conditions 
for the airfoil model and the upper and lower walls are independent 
of z. It is also assumed that the tunnel is narrow enough for the flow 
at each sidewall to be strongly Influenced by the other sidewall bound- 
ary layer. Reference 6 indicates that to the lowest order, the spanwlse 
velocity varies linearly with the spanwlse coordinate z as 


_ 2Uz 36* 
b 9x 


(3) 


In wider tunnels the disturbance caused by the sidewall boundary layer 
decays nonlinearly with distance from the sidewall so that equation (3) 
is not valid. 

The flow in the wind tunnel described above is governed by the 
small perturbation form of the continuity equation, which can be written 
as 


(l . + + 

' «> ' 3x dy oz 


(Y H |1 |!1 

“ U 3x 


( 4 ) 
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where is the freestream Mach number, u ■* U - is the velocity 

perturbation in the x-direction, and y is the ratio of specific heats 
of the gas. 

The dynamics of the sidewall boundary layer are modeled with the 
von Karman momentum integral, which is given in reference 12 and can be 
written as 


96 * 

3x 


-^(2 

U ' 


+ u ^ 


M 4. _!». 

pU^ 


(5) 


where p is the density and 6 , and H are the sidewall dis- 

placement thickness, wall shear stress, and shape facto:, respectively. 

For the present problem, equation (5) can be simplified because the 
sidewall boundary layers in the test sections of most airfoil wind 
tunnels can be approximated as flat-plate boundary layers with large 
equivalent lengths, L, and hence, relatively large Reynolds nximbers. 

The model pressure field is considered to cause a rather localized 
variation in this large length-scale sidewall boundary layer, and by 
applying the following order of magnitude analysis, the shear stress 
term can be neglected from equation (S). 

3H 

First, the shape factor gradient, can be related to the 

3u 

velocity gradient, by the following expression derived in 

Appendix B: 


. (H 4- 1)(H - 1) ^ 
9x U 9x 


( 6 ) 
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Therefore, for a flat plate boundary layer without a pressure gradient, 

both 3U/3x and 3H/3x vanish in equation C5), leaving only the shear 

stress tern, T^/pU^, to affect the sidewall boundaty-layer growth rats, 

35*/3x. The order of nagnitude of 35*/3x in the test section is 

which should be the same order as T«/pU^. 

For the sidewall boundary layer with a pressure gradient due to 
dU 

the model, — is of the order ^ so that the first two terms in 
equation (3) are of the order 5*/c. Because the equivalent length of 
the sidewall boundary layer, L, is much larger than the model chord 
length, c, the inequality 




(7) 


j. 2 

applies^, and the shear stress term, neglected from 
equation (3) as a first approximation. With equation C6), the final 
form becomes 




(2 + 


1 - M 

H ^ 3x 


( 8 ) 


With equation C8), equation (3), and the observation that 


3u 

3x " 3x 


the derivative 




( 9 ) 


'Measurements of skin friction and 6 showed values of T^/pu 
ranging from 0.0010 to 0.0012 and values of 5*/c ranging from 0.014 
to 0.052, which experimentally verified inequality (7). 
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is obtained. For small disturbance flow, H in equation (9) can be 
replaced by the freestream Hach number Equation (9) can be com- 
bined with equation (4) to give 




(Y + 


U ^ 


( 10 ) 


or 


^ 3y 


(Y 


^ u ^ 


where 
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For this study, the values of 5 and U measured at the model 
station, tunnel empty, are used as constant values in equation (10) so 
that the von Karraan transonic similarity rul^, discussed in reference 15, 
can be applied. This rule relates the pressure coefficients of two flow 
fields, denoted by Cp and Cp 2 < 


<T ^ D<1 


(Y + 1)M„ 


2 


'Ptl 
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'P.2 


(13) 


where the flow fields satisfy the constraint 
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For the same model and test gas (t^ - tj* “ Y 2 )* equation (14) 
becomes 
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6 . 


6 . 


(15) 


An interference-free equivalent Mach number can be defined 

it 

with equation (IS) and the condition 5 « 0 as 



M 


00 



(16) 


This equivalent Mach number represents a flow In a: ideal two- 
dimensional tunnel without a sidewall boundary layer which is otherwise 
the same as the actual two-dimensional tunnel with a sidewall boundary 
layer. The pressure coefficient can be adjusted from the value in the 
actual two-dimensional tunnel to the value in the ideal two-dimensional 
tunnel having a similar flow without a sidewall boundary layer with 
equation (13) with Yj “ Y 2 (^aoe test gas). The expn.jsion is 




(17) 


where Cp is the pressure coefficient in the ideal two-dimensional 
tunnel. Equation (17) results in the following adjusted airfoil force 
coefficients formulated by integrating the surface pressures: 
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adjusted section normal-force coefficient, c^. 



adjusted section drag coefficient, c^. 


IS 


(18) 



where c_ and 
n 

coefficients. 


6 



c^ are the measured section normal-force and drag 


(19) 


Approximate Mach Number Increment 
An approximate expression for the increment can be 

formulaced from a firsc-order Taylor series expansion of both sides of 

m • 

equation (16). First, and B are rewritten as. 


+ AM 

OO 00 


( 20 ) 


and 


/Tv 


+ A6 


( 21 ) 


where 


AB 




( 22 ) 


and is treated as a single variable. Equation (16) then appears as 


M + AM 

GO 



(M^ + AM) 




3/ii 


M 


.1 - y 


.3/A 


(23) 


I 


1 



f 


ORIGINAL PAGE 18 
or POOR QUALITY 


where the only variables are AM and A6> Both sides are expanded 
about AM and A3 equal to sero. Retention of the first-order terns 
of both series results in 


M_ + AM 


1 , 3 “ 




771 * 


5 M 


3 

374 ■ 4 


AB 


(1 - (• - 


7/4 


(24) 


The first-order approxintatlon is 

3M AB 
00 

4(1 - 

\ CO 

This equation can be solved for AM in terms of A6 to give 




(25) 


AM 



(26) 


or 


AM 




3M 2 + 




2 

2 + M 

00 


(27) 


For M^ ranging between 0.7 and 0.9 and H ranging from 1.4 to 1.6, 
this Increment is approximately 



which represents the fraction of the tunnel width occupied by the two 
sidewall boundary-layer displacement thicknesses. 
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CHAPTER 4 

EXPERIMENTAL RESULTS 
Equivalent Freestream Mach Nuaber 

The effectiveness of the sidewall similarity rule has been evaluated 
by comparing measured airfoil test data obtained at different sidewall 
boundary-layer displacement thicknesses. First, the equivalent Mach 
number, 3^, was determined for the three boundary-l?yer thicknesses 
with equations (12) and (16). The value- of S* and H used in equa- 
tion (12) were measured at the model station, tunnel empty, as suggested 
by the analysis. Figure 10 shows the variation of 3^ with for 

the three boundary-layer thicknesses and shows an increment between the 
equivalent Mach number and the measured freestream Mach number of 
approximately 26*/b as Indicated in equation (28). 

NASA 0012 and SC-27 Airfoil Tests 

Next, the airfoil tests were conducted beginning with the NACA 0012 
airfoil. This airfoil has some interesting transonic behavior at zero 
angle of attack in that the chordwise shock wave location varies almost 
linearly with freestream Mach number up to values of approximately 0.84. 
The variations in shock wave location wit.i both and 3^ were com- 

pared for the three sidewall boundary-layer thicknesses, as shown in 
figure 11. A significantly improved correlation was obtained when 
rather than N was used. 

CO 

The data for the first two sidewall boundary- layer thicknesses are 
shown in figure 11 with centered symbols, while the data for the third. 
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or thickest* sidewall boundary layer are shown with open symbols. Other 
observations that will be awntioned later in the paper indicated that th^ 
thickest sldemll boundary layer aay have introduced excessive three 
diaensionallty that was not addressed by the present analysis. There- 
fore, the centered symbols denote the data for which the similarity rule 
can be applied with the most confidence. 

The other two transonic characteristics investigated were the 
variation of normal-force and section drag coefficient with freestream 
Mach number at a fixed angle of attack. Here, the similarity rule 
requires the application of equations (18) and (19) to form the values 
of the adjusted normal-force and section drag coefficients* c^ and c^. 
The variations of the measured normal-force and section drag coeffi- 
cients, c_ and c., with the measured freestream Mach number, M , 
no «> 

were compared to the variations of the adjusted normal-force and section 
drag coefficients, c^ and c^, with the equivalent freestream Mach 
number, M^, for the three sidewall boundary- l"yer thicknesses. 

Continuing with the NACA 0012 airfoil at zero angle of attack, this 
phase of the investigation began with the comparison of the variation of 
the measured drag coefficient, c^, wi^h M^ to the variation of the 
adjusted drag coefficient, c^, with M^. This adjustment actually 
applies only to the component of pressure drag in the drag coefficient 
and does not account for the skin-friction component. Figure 12 shows 
the comparison between c^ vs. M^ and Cj vs. for the three side- 
wall boundary layers. In figure 12, the similarity rule provides a 
substantially improved drag correlation in the drag-rise region, but 
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loses qtiallty below the drag rise. This is probably because the 
majority of the drag comes from the skin friction below drag rise 
whereas the adjusted drag coefficient addresses only the pressure drag. 
The correlation improves as the pressire drag becomes a larger fracr^rn 
of the total drag, as seen in the drag rise region. 

Figure 12 indicates more scatter in the drag data measured with tl» 
thickest sidewall boundary layer. This boundary layer was approximately 
5.2 cm thick at the model station, tunnel enpty, so that che two side- 
wall boundary layers occupied approximately two-thirds of the tunnel 
width. This large amount of sidewall boundary layer is thought to cause 
three-dimensional secondary flows not addressed by the analysis, and 
could possibly have adverse effects on both the drag measuren'ents with 
che wake probe and the airfoil su.face pressure measurements. There- 
fore, the data for the thickest sidewall boundary layer are presented 
with open symbols, while the data for the two thinner sidewall boundary 
layers, where Che similarity rule is more applicable, are presented with 
centered symbols. This depiction was also used in figure 11 and will 
follow for all figures presenting measured airfoil data with the three 
sidewall boundary- layer thicknesses. 

The variation of normal-force coefficient with freestream Mach num- 
ber was first studied using the NACA 0012 airfoil at one degree angle of 
attack. Again, the similarity rule was evaluated by comparing Che 
variation of the measured normal-force coefficient, c^, with to 

the variation of the adjusted normal-force coefficient, c^, with 1^. 
Figure 13 shows this comparison. The similarity rule provides a 
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significantly improved correlation, particularly for the two thinner 
sidewall boundary layers. The correlation quality diminishes for the 
thickest sidewall boundary layer, probably because of the previously 
mentioned problons associated with its large thickness. Some loss in 
correlation also appears for the two thinner sidewall boundary layers 
near the maximum normal-force coefficient. This is thought to be caused 
by interactions between the shock wave on the airfoil upper surface with 
the sidewall boundary layer. 

The variation of section drag coefficient with freestream Mach num- 
ber was also studied for the NACA 0012 airfoil at one degree angle of 
attacK. Figure 14 provides a comparison of vs. and c^ vs. 3^. 
The results are generally similar to the zero angle-of -attack case for 
this same airfoil in that the correlation improves in the drag rise 
region. 

The SC-27 supercritical airfoil was used to examine the correlation 
of the adjusted normal-force coefficient with the equivalent Mach number. 
This airfoil has a much weaker shock wave than that on the NACA 0012 air- 
foil at the same normal-force coefficient. This would hopefully alle- 
viate the possible interaction between the airfoil shock wave and the 
sidewall boundary layer that was suspected for the NACA airfoil. 

Figure 15 shows the comparison between c^j vs. and c^ vs. 
for the SC-27 airfoil at zero angle of attack. The correlation is very 
good for the first two sidewall boundary layers, but degraded somewhat 
for the thickest sidewall boundary layer. In comparison with the 
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NACA 0012 airfoil shown in figure 13, the correlation for the SC-27 air- 
foil for the two thinnest sidewall boundary layers is improved even with 
the SC-27 at higher values of adjusted normal-force coefficient. 

The correlation of c with M was examined for the SC-27 airfoil 

n 00 

at higher values of c^ and, therefore, stronger airfoil shock waves. 

Figure 16 shows the comparison between c^ vs. and c^ vs. for 

an angle of attack of one degree. While the c^ vs. correlation is 

still much better than that for c vs. M , the correlation for the two 

n oo» 

thinnest sidewall boundary layers was of lower quality than that shown 
in the zero angle-of-attack case in figure 15. Therefore, the shock 
wave strength on the airfoil appears to limit the performance of the 
similarity rules. 


NLR-1 Airfoil Test Results 

The NLR-1 airfoil test provided a study of the correlation of the 
adjusted normal-force and section drag coefficients with equivalent Mach 
numbers for several angles of attack. As mentioned in the models' 
description, the NLR-1 is a supercritical airfoil for rotorcraft 
applications. 

Figure 17 shows the comparison of c^ vs. and c^j vs. for 
an angle of attack of zero. The correlation for the three sidewall 
boundary layers is improved using c^ vs. rather than c^ vs. 

but with different results than those observed for the NACA 0012 and 

m 

SC-27 airfoil data. The largest values of for the second sidewall 

boundary layer exceed those of the first or thinnest sidewall boundary 


layer. 



1 


22 


The reasoa the second sidewall boundary layer produced a larger 
maximum value of c^ than the first sidewall boundary layer was deter- 
mined from an examination of the local Mach number distributions. 

Figure 18(a) shows the local Mach number distributions on the hXR-1 
airfoil with the two thinnest sidewall boundary layers. The adjusted 
normal force coefficients are near their mavlmiim respective values ^own 
in figure 17, and these values occur at essentially the sane equivalent 
freestream Mach number. Except for the lower-surface region near 
30 percent chord, the local Mach numbers for the thinner sidewall 
boundary layer are slightly less than those for the thicker sidewall 
boundary layer at the same location. This condition is required for 
matched values of the adjusted pressure coefficient which relate 

directly to the value of c . The value of c for the thinner side- 

n n 

wall boundary layer is slightly lower than that if the thicker sidewall 
boundary layer because the lower-surface local Mach numbers at 30- 
percent chord for both sidewall boundary layers are practically the sane. 
This causes the adjusted pressure coefficients for the thinner sidewall 
boundary layer to have a larger negative magnitude than that of the 

thicker sidewall boundary layer in this lower-surface region, and 

■ 

results in a lower value of c . 

n 

Figure 18(b) compares the local Mach number distributions for the 
same two sidewall boundary layers with a small increase in equivalent 
freestream Mach number. With the thinner sidewall boundary layer, the 
local Mach numbers on the lover surface between 20- and 50-percent chord 
have substantially increased from those seen in figure 18(a) with only a 
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very small change in the equivalent freeatream Mach number (0.850 to 
0.864). These local Mach numbers for the thinner sidewall boundary 
layer are significantly higher than those for the thicker sidewall 
boundary layer. The corresponding higher negative magnitude of the 
adjusted pressure coefficients in this lower-surface region produce a 
much lower value of c^ for the thinner sidewall boundary layer. The 
local Mach numbers for the thicker sldemll boundary layer in this same 
lower-surface region also have a noticeable increase in values compared 
to those in figure 18(a). The equivalent freestream Mach number for 
the thicker sidewall boundary layer has only changed from 0.852 
(figure 18(a)) to 0.858 (figure 18(b)). This sensitive development of 
supersonic flow on the lower surface caused an abrupt loss in normal 
force with Increased equivalent freestream Mach number. 

This study included two positive angles of attack shown in figures 
21 and 22. The correlations of the normal-force coefficient with 
freestream Mach number for these angles of attack appear very similar 
to the results for the NACA 0012 airfoil at one degree angle of attack, 

S ■ 

r.s indicated by figure 13. First, the vs M^ provides an 

unquestionable Improvement in correlation over the c^ vs M^. Second, 
the correlation of c vs M loses quality for the data with the 

n 00 

thickest sidewall boundary layer as compared to the two thinner side- 
wall boundary layers. Third, these two thinner sidewall boundary layers 
show a slight loss in correlation at the maximum values of c^, which 
is probably due to the presence oi strong shock waves on the airfoil 
interacting with the sidewall boundary layer. These shock waves 
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contributed largely to the rapid drop In that follows the 

naxlBum value. 

Drag aeasurensnts were also obtained for the NLR-1 airfoil. These 
drag data, which are presented in figures 23 through 26, in nany 
cases indicate very similar correlation behavior to that shown for 
the NAGA 0012 airfoil In figures 12 and 13. The correlation of 

m m 

Cj vs shoirs no real Improvemoit over c^ vs until the drag 

rise region. Drag measurements Involving the thickest sidewalJl 
boundary layer often seem to show more scatter than that from the 
t%ro thinner sidewall boundary layers, again, probably because of the 
large fraction of the tunnel width occupied by the thickest sidewall 


boundary layer. 
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Sunaary of Airfoil Tests Results 

For the airfoil data measured with the three ^tidewall boundary 
layers, the similarity rule produced an improved correlation for the 
variation of the adjusted normal-force and section drag coefficients 
with the equivalent freestream Mach number as compared to the variation 
of the measured coefficients with the measured freestream Mach number. 

The normal-force coefficients appear to form three distinct zones 
for each sidewall boundary layer when plotted against the measured free- 
stream Mach number. The adjusted normal-force coefficients appear to 
have more converged zones for the three sidewall boundary layers when 
plotted against the equivalent freestream Mach number, but the magni- 
tudes of the adjusted normal-force coefficients do not entirely 
converge . 

The section drag coefficients appear to have two distinct drag rise 
regions for the two thinnest sidewall boundary layers when plotted 
against the measured freestream Mach number. The adjusted section drag 
coefficients show converged drag rise regions when plotted against the 
equivalent freestream Mach number but show no improvement in the magni- 
tude of the drag coefficient when below drag rise. 

An important effect of the similarity rule is that the maximum 
adjusted normal-force coefficient and the divergence that follows occur 
at almost the same equivalent freestream Mach numbers for all three 
sidewall boundary layers. Likewise, the drag rise occurs at almost the 
same equivalent Mach number for the two thinnest sidewall boundary 
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layers. The fact that the maxlmuffl nomal force and drag rise occur 
at different measured freestream Mach ntanbers fcr each of the three 
sidewall boundary layers but at the same equivalent freestream Mach 
number demonstrates the correction to the measured freestream Mach 
number provided by the similarity rule. 
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CHAPTER 5 

CONCLUDING REMARKS 

The effects of attached sidewall boundary layers In two-dimensional 
transonic tunnels have been correlated with a transonic similarity rule. 
It ha been shown experimentally that the application of this similarity 
rule to the airfoil test data obtained in the Langley 6- by 19-Inch 
Transonic Tunnel gives an effective freestream Mach number correction. 
The experimental data also Indicate that the similarity rule provides a 
substantial correction to the normal-force coefficients and some correc- 
tion for the section drag coefficients in the drag rise region. 

The similarity rule correction applies provided the sidewall 
boundary layer is small enough tu avoid excessive three-dimensional 
interactions with the model. The similarity rule can be used as long as 
the sidewall boundary layers have no appreciable separation (due to 
shock wave/boundary layer interaction or significant trailing-edge 
separation) . 
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APPENDIX A 

BOUNDARY-LAYER DATA REDUCTION 
Temperature Distribution in the Boumlary Layer 


The velocities within the boundary layer were calculated from the 
values of local Mach number and speed of sound. The local temperature 
within the boundary layer was given by equation (1)» where the wall 
temperature was assumed to be the adiabatic temperature indicated in 
reference 13 as 


aw 


T 1 + r 




(Al) 


For the tunnel-empty sidewall boundary layer it was assumed that 
- Tg, « Mg, and the recovery factor, r, was given a value of 
0.8963, which was also obtained from reference 13. The temperature of 
the chamber surrounding the test section usually remained within 5°K of 
the value for in equation (Al). Due to the rapid operation in 

blowdown testing, the wall temperature was not expected to change 
signif icantly. 

The temperature distribution through the boundary layer was obtained 
from reference 12 as 


T “ T + (T 

aw 'aw 


- T ) i 

‘w^ Ug 



(A2) 



where f Is the time-averaged local temperature at some point within 

the boundary layer. Since is assumed to be the actual wall 

2 2 

tmsperature, and ■ M^Y(Cp ~ c^)Tg, equations (Al) and (A2) can 

be ccHDbined to provide equation (1) in the text. 



where T is assumed to be the local static temperature, T. 


(A3) 


(A4) 


Displacement Thickness and Momentum Thickness Calculations 
The sidewall boundary- layer thickness was determined from the 
least-squares power-law regression given by equation (2) in the text. 
This procedure was used because the largest sidewall boundary- layer 
thickness in the experiment exceeded the highest total head tube on the 
boundary-layer rake-tube probe. Using the power-law representation of 
the velocity profile in the boundary layer allowed a simple, closed-form 
integration for calculating the displacement thickness and the momentum 
thickness. 

The boundary-layer displacement thickness 6 is defined as 
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From the ideal gas relation, p ° pRT, and the assumption of zero normal 

p "^e 

pressure gradient in the boundary layer, the expression ~:jr is 

Pe 1 

obtained so that equation (A4) becomes 


_e_ 


1 + 0.1793 


4-S 


(A6) 


Using the power-law representation of the velocity profile given in 
equation (2) and equation (A6), equation (AS) becomes 


5* = 


f 


1 - 


(z/c) 


1/N 


dz 


1 + 0.1793 M^[l - (z/5)^^^ 

Note that the second term of the integrand is of the form 


(A7) 



~ 1 ” 

1 + 0.1793 

e 

, , ,.,2/N 

1 - o(z/6) 

1 + 0.1793 M 

e _ 


(A8) 


which can be rewritten as the geometric series 


— 5-[i + P(z/6)^^*^ + P^(z/6)^^” + P^(z/6)^^” ..] 

1 + 0.1793 M 

e 

(A9) 

where 


0.1793 

P — ^ (AlO) 

1 + 0.1793 M 


1 
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Placonent of the series in expression (A9) into equation (A7) allows a 
term by term integration which results in the series 


6 * - 6 


It ^ ( 

N NP Np2 \ 

^ ^ 2\I 

1 + 0.1793 MV 
_ e 

4 + 1N + 3 N + 5 ‘ ) 


(All) 


or 


6 * = 6 + 


1 + 0.1793 M 


“ k-1 

NP* ^ 


X) (2k - 1) + N 
k=l 


(A12) 


This series converges rapidly and yields the necessary precision when 

k - 5. 

The momentum thickness is defined as 



(A13) 


and is calculated in a manner similar to that used for the displacement 
thickness. With the power-law representation of the velocity profile 
and the definition of 6*, equation (A13) becomes 


0 “ 6 - 



(A14) 


The last term in this equation can also be represented by the geometric 
series 


. pU/5)2« 

1 + 0.1793 

6 




■] 


(A15) 
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0.1793 M 


irtiere P » 


1 f 0.1793 


as in equation (AlO) . Term by term 


integration yields 


6 r N NP NP^ NP^ “I 

1 + 0.1793 + 2 N + 4 N + 6 N + 8 -J 


1 + 0.1793 M‘ 


“* k-1 

E NP 

2k + N 


(A16) 


which results in the final equation 


9 - 6 - 5 * - 


1 + 0.1793 M 


NP 


k-1 


2k + N 


e k=l 


(A17) 


This summation also required 5 terms for convergence to the necessary 


precision. 
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APPENDIX B 

THE RELATIONSHIP OF THE VELOCITY GRADIENT TO 
THE SHAPE FACTOR GRADIENT 

The conventional shape factor, H, has a significant dependence on 
freestream Mach number and, therefore, is often replaced by the trans- 
formed shape factor, H. Reference 16 defines the transformed shape 
factor as 



(Bl) 


For compressible, turbulent boundary layers with constant total tempera- 
ture assumed through the boundary layer, reference 16 indicates that H 
is related to H by 

h = (H + 1)^1 + ^ 

Reference 16 also shows that for large Reynolds numbers, such as Chose 

_ 

applicable to tunnel sidewall boundary layers, H approaches one. 

This simplifies equation (B2) to 

H - 1 + (Y - 1)M^ (B3) 


The measured values of H ranged from 1.18 to 1.26, but use of 
these values in the above analysis did not provide any significant 
difference from using H - 1. 
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Use of the simple compressible flow relations with coni-tant total 
temperature results in the expression. 


~ m\rt 

.2 o 


1 + 




(B4) 


Differentiation of both sides with respect to x gives 


2u|!^ - 

9x 


2MYRT 


M YRI^Cy - 1)M 





(B5) 


Division of both sides by 2U‘^ and use of equation (B4) yields 


1 ^ 
U 9x 


.21 / 


m[ 2 + (Y - 1)M^J 
which relates the velocity gradient to the Mach number gradient. 

The shape factor gradient can be related to the Mach number 
gradient by 


3H _ 9M 

^ - 20 - DM ^ 


(B6) 


(B7) 


with equation (B3), and the Mach number gradient can be related by the 
velocity gradient with equation (B6), so that 


|S- (Y - Um2[2^ (Y - 


(B8) 


2 

With equation (BJ), the expression (y - 1)M - H - 1 results, and the 

final form of equation (6) in the text is obtained as 


9H _ (H -b 1)(H - 1) 9U 
9x U ?x 


(B9) 
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TOP VIEW 
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Figure 4,- Experimental apparatus used in the Langley 6- by 19-Inch 
Transonic Tunnel to investigate the effects of the sidewall 
boundary-layer displacement thickness on two-dimensional testing. 
All dimensions are in cm. 
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(a) Top and front view. 

Figure 5.- Sketch of the rake-tube probes used to survey the sidewall 
bounaary layers in the Langley 6- by 19-Inch Transonic Tunnel 
All dimensions are in cm. 
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Figure 6.- Nondimens ional velocity distribution in artificially 
thickened sidewall boundary layers. 
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A.O.A.= 0.0 deg.: FIXED TRANS, oi 0.075 x/c; RN=3.4-3.8 X 10* 
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(a) Shock wave location vs. measured freestreeun Mach number. 

Figure 11.- Variation of shock wave location with freestream Hach 
number for the NACA 0012 airfoil tested with three sidewall 
boundary- layer displacement thicknesses. Angle o£ attack is 
0 degrees. 
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A.O.A.= 0.0 deg.: FIXED TRANS, at 0.075 x/c; RN=3.4-3.8 X 10* 
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(a) Section drag coefficient vs. measured freestream Mach number. 

Figure 12.- Variation of section drag coefficient with freestream Mach 
number for the NACA 0012 airfoil tested with three sidewall boundary- 
layer displacement thicknesses. Angle of attack is I degrees. 
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A.O.A.= 0.0 deg.; FIXED TRANS, at 0.075 x/c; RN=3.4-3.8 X 10® 
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(b) Adjusted section drag coefficient vs. equivalent 
freestream Mach number. 

Figure 12.- Concluded. 
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(a) Normal-force coefficient vs. measured freestream Mach number. 

Figure 13.- Variation of normal-force coefficient with freestreeun 
Mach number for the NACA 0012 airfoil tested with three siCewall 
boundary- layer displacement thicknesses. Angle of attack is 
1.0 degree. 
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A.O.A.= 1.0 deg.; FIXED TRANS. o1 0.075 x/c; RN=3.4-3.8 X 10* 
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(b) Adjusted normal-force coefficient vs. equivalent freestream 

Mach number. 


Figure 13.- Concluded. 
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A.O.A.= 1.0 deg.; FIXED TRANS, ot 0.075 x/c; RN=3.4-3.8 X 10* 
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(a.) Section drag coefficient vs. measured freestreeun Mach number. 

Figure 14.- Variation of section drag coefficient with freestream 
Mach number for the NACA 0012 airfoil tested with three sidewall 
boundary- layer displacement thicknesses. Angle of attack is 
1.0 degree. 
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A.O.A.= 1.0 deg.; FIXED TRANS, at 0.075 x/c; RN=3.4-3.8 X 10* 
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(b) Adjusted section drag coefficient vs. equivalent freestream 

Mach number. 


Figure 14.- Concluded. 
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A.O.A.= 0.0 deg.; FIXED TRANS, ot 0.050 x/c; RN=5.4-3.8 X 10* 
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(a) Normal-force coefficient vs. measured freestreeun Mach number. 

Figure 15.- Variation of normal-force coefficient with freestream 
Mach number for the SC-27 airfoil tested with three sidewall 
boundary- layer displacement thicknesses. Angle of attack is 
0 degrees. 
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(b) Adjusted normal-force 'efficient with equivalent freestream 

number. 

Figure 15.- Concluded. 
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A.O.A.= 1.0 deg.; FIXED TRANS, a* 0.050 x/c; RN=3.4-3.8 X 10* 
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(a) Normal-force coefficient vs. measured freestream Mach number. 

Figure lb.- Variation of normal-force coefficient with freestreeun 
Mach number for the SC-27 airfoil tested with three sidewall 
boundary- layer displacement thicknesses. Angle of attack is 

1.0 degree. 
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A.O.A.= 1.0 deg.; FIXED TRANS. ol 0.050 x/c; RN=3.4-3.8 X 10* 
2d’/b 





(h' Adjusted nonaal-force coefficient vs. equivalent freestream 

Mach number. 




Figure 16.- Concluded. 
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(a) Nonoal'force coefficient vs. measured freestream Mach number. 

Figure 17.- Variation of normal-force coefficient with freestream 
Mach number for the NLR-1 airfoil tested with three sidewall 
boundary- layer displacement thicknesses. Angle of attack is 
0 degrees. 
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(b) Adjusted normal-force coefficient vs. equivalent freestream 

Mach number. 


Figure 17.- Concluded. 
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LOCAL MACH NUMBER DISTRIBUTIONS ON NLR-1 AIRFOIL 
A.O.A^ 0.0 DEG.. FIXED TRANS. AT1o.05 X/G 
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Figure 18.- Chordvlse local Hach number distributions on the 
NLR-1 airfoil. Angle of attack is 0.0 degrees. Open symbols 
Indicate the airfoil upper surface; centered symbols Indicate 
the airfoil lower surface. 
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A.O.A.= 0.0 DEG., FIXED TRANS. AT 0.05 X/C 

o M^rO.864. 2«Vb =0.028. RN=3.64 X 10*. Cp = 0.045 
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(b) 5^ = 0.86 

Figure 18.- Concluded. 
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A.O.A.=-1.0 deg.; FIXED TRANS, oi 0.050 x/c; RN=3.4-3.8 X 10* 
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(a) Normal-force coefficient vs. measured freestream Mach number. 

Figure 19.- Variation of normal-force coefficient with freestream 
Mach number for the NLR-1 airfoil tested with three sidewall 
boundary- layer displacement thicknesses. Angle of attack is 
-1.0 degree. 
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(b) Adjusted normal-force coefficient vs. equivalent freestream 

Mach number. 


Figure 19.- Concluded. 
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LOCAL MACH NUMBER DISTRIBUTIONS ON NLR-1 AIRFOIL 
A.O.A.=-1.0 DEG., FIXED TRANS. AT 0.05 X/C 



O m^= 0.853. =0.028. RN=3-72 X 10*. Cp = -0.044 

o m^= 0.858. 2a*/b =0.070, RN=3.90 X 10*. ^^=-0.003 
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Figure 20.- Chordwlse local Mach number distributions on the 
NLR-1 airfoil. Angle of attack Is -1.0 degree. Open svnbols 
Indicate the airfoil upper surface; centered symbols inolcate 
the airfoil lower surface. 
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A.O.A.= 1.0 deg.; FIXED TRANS, at 0.050 x/c; RN=3. 4-3.8 X 10* 
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(a) Normal-force coefficient vs. measured freestream Mach number. 

Figure 21.- Variation of normal-force coefficient with freestream 
Mach number for the NLR-1 airfoil tested with three sidewall 
boundary- layer displacement thicknesses. Angle of attack is 
1.0 degree. 
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A.O.A.= 1.0 deg.; FIXED TRANS, ai 0.050 x/c; RN=3.4-3.8 X 10* 
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(b) Adjusted normal-force coefficient vs. equivalent freestream 

Mach number. 


Figure 21.- Concluded. 
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A.O.A.= 2.0 deg.; MXED TRANS, ai 0.050 x/c; RN=3.4-3.8 X 10* 
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(a) Normal-force coefficient vs. measured freestream Mach number. 

Figure 22.- Variation of normal-force coefficient with freeitream 
Mach numoer for the NLR-1 airfoil tested with three sidev. 11 
boundary- layer displacement thicknesses. Angle of attack is 
2.0 degrees. 
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(b) Adjasted normal-force coefficient vs. equivalent iresstream 

Mach number. 


Figure 22.- Concluded. 
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A.O.A.= 0.0 deg.; FIXED TR^;iS. at 0.050 x/c; RN=3.4-3.8 X 10* 
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(a) Section drag coefficient vs. measured freestream Mach nuinber.. 

Figure 23.- Variat.ion of section drag coefficient with freestream 
Mach number . tne NLR-1 airfoil tested with three si'^'^'/all 
boundary- lay .r displacement thicknesses. Angle of attack is 
0 degrees. 
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A.O.A.= 0.0 deg.; FIXED TRANS, at 0.050 x/c; RN=3.4>3.8 X 10* 
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(b; Adjusted section drag coefficient vs. equivalent freestream 

Mach number. 


Figure 23.- Concluded. 
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A.O.A.=-VO deg.; FIXED TRANS, at 0.050 x/c; RN=3.4-3.8 X 10* 
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(a) Section drag coefficient vs. measured freestream Mach number. 

Figure Variation of section drag coefficient with freestrofun 

Mach number for the NLR-1 airfoii. tested with three sidewall 
boundary- layer displacement thicknesses. Angle of attack is 

-1.0 degree. 
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(bj Adjusted section drag coefficient vs. equivalent freestream 

Mach number. 




Figure 24.- Concluded. 
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(a) Section drag coefficient vs. measured freestream Mach number. 

Figure 25 .- Variation of section drag coefficient with freestream 
Mach number for the NLR-1 airfoil tested with three sidewall 
boundary- layer displacement thicknesses. Angle of attack is 
1.0 degree. 
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(b) Adjusted section drag coef .'icient vs. equivalent freestream 

Mach nuuber. 


Figure 25.- Corcluded. 
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A.O.A.= 2.0 deg.; FIXED TRANS, at 0.050 x/c; RN=3.4-3.8 X 10* 



(a) Section drag coefficient vs. measured freestream Mach number. 

Figure 26.- Variation of section drag coefficient with freestream 
Mach number for the NLR-1 airfoil tested with three sidewall 
boundary- layer displacement thicknesses. Angle of attack is 
. 0 degrees . 
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(b) Adjusted section drag coefficient vs. equivalent freestream 

Mach number. 

Figure 26.- Concluded. 
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